Метод степеневого ряду використовується для пошуку розв'язку у вигляді степеневого ряду диференціального рівняння. В загальному, цей підхід розглядає степеневий ряд з невідомими коефіцієнтами і підставляє його в диференціальне рівняння, щоб знайти рекурентне співвідношення для коефіцієнтів.
Розглянемо лінійне диференціальне рівняння другого порядку
![{\displaystyle a_{2}(z)f''(z)+a_{1}(z)f'(z)+a_{0}(z)f(z)=0.\;\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e36950769e73bec96b189c3f16b05de752760bb1)
Припустимо, що a2 не нуль для всіх z. Тоді ми можемо поділити і отримати
![{\displaystyle f''+{a_{1}(z) \over a_{2}(z)}f'+{a_{0}(z) \over a_{2}(z)}f=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/51d51f25f5571ee6df1d2b463a6a0957d88ecdc0)
Припустимо також, що a1/a2 і a0/a2 — це аналітичні функції.
Метод степеневого ряду потребує побудови степеневого ряду розв'язку
![{\displaystyle f=\sum _{k=0}^{\infty }A_{k}z^{k}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f37bf5646d2e28f8eac837567f7c5f0d6c17377f)
Якщо a2 обертається в нуль для деякого z, тоді можна викоритсати метод Фробеніуса, видозміна цього методу, для знаходження розв'язку поблизу сингулярностей. Метод також працює для рівнянь більш високого порядку або систем.
Давайте розглянемо Ермітове диференціальне рівняння,
![{\displaystyle f''-2zf'+\lambda f=0;\;\lambda =1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f466a38f18c086d67c94ea660a391c0e368a8b8)
Ми можемо спробувати побудувати ряд розв'язку
![{\displaystyle f=\sum _{k=0}^{\infty }A_{k}z^{k}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d40af09eea80208a2a1c77178387b7dca981224)
![{\displaystyle f'=\sum _{k=1}^{\infty }kA_{k}z^{k-1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7647fc96eb617f9605e816a60f56271e39c8b5b)
![{\displaystyle f''=\sum _{k=2}^{\infty }k(k-1)A_{k}z^{k-2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5bf1d5701d5a83dbc2bdf84d38890aecd9fe1730)
Підставляючи це в диференціальне рівняння
![{\displaystyle {\begin{aligned}&{}\quad \sum _{k=2}^{\infty }k(k-1)A_{k}z^{k-2}-2z\sum _{k=1}^{\infty }kA_{k}z^{k-1}+\sum _{k=0}^{\infty }A_{k}z^{k}=0\\&=\sum _{k=2}^{\infty }k(k-1)A_{k}z^{k-2}-\sum _{k=1}^{\infty }2kA_{k}z^{k}+\sum _{k=0}^{\infty }A_{k}z^{k}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37f8afe1b10108e42d6197bd1d30af30dbfbe10d)
Зсуваючи перший доданок
![{\displaystyle {\begin{aligned}&=\sum _{k=0}^{\infty }(k+2)(k+1)A_{k+2}z^{k}-\sum _{k=1}^{\infty }2kA_{k}z^{k}+\sum _{k=0}^{\infty }A_{k}z^{k}\\&=2A_{2}+\sum _{k=1}^{\infty }(k+2)(k+1)A_{k+2}z^{k}-\sum _{k=1}^{\infty }2kA_{k}z^{k}+A_{0}+\sum _{k=1}^{\infty }A_{k}z^{k}\\&=2A_{2}+A_{0}+\sum _{k=1}^{\infty }\left((k+2)(k+1)A_{k+2}+(-2k+1)A_{k}\right)z^{k}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f444f79db3a0bc28d5f710967268abac65a15382)
Якщо цей ряд є розв'язком, тоді всі ці коецфіцієнти дорівнюють нулю, отже для k=0 і k>0:
![{\displaystyle (k+2)(k+1)A_{k+2}+(-2k+1)A_{k}=0\;\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0ac743f444065cf020f3e80abc320c6d2877ebd9)
З цього ми можемо отримати рекурентне співвідношення для Ak+2.
![{\displaystyle (k+2)(k+1)A_{k+2}=-(-2k+1)A_{k}\;\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/39c5487720c093f6a146062740bced2830e54cd1)
![{\displaystyle A_{k+2}={(2k-1) \over (k+2)(k+1)}A_{k}\;\!}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc85f067565b87d1c74c285b2d1e3083b4e4cda3)
Тепер маємо
![{\displaystyle A_{2}={-1 \over (2)(1)}A_{0}={-1 \over 2}A_{0},\,A_{3}={1 \over (3)(2)}A_{1}={1 \over 6}A_{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/41745fa16634146049088e03e633e575e91a9100)
Ми можемо знайти A0 і A1 якщо нам задані початкові умови, тобто якщо ми маємо задачу Коші.
Отже, маємо
![{\displaystyle {\begin{aligned}A_{4}&={1 \over 4}A_{2}=\left({1 \over 4}\right)\left({-1 \over 2}\right)A_{0}={-1 \over 8}A_{0}\\[8pt]A_{5}&={1 \over 4}A_{3}=\left({1 \over 4}\right)\left({1 \over 6}\right)A_{1}={1 \over 24}A_{1}\\[8pt]A_{6}&={7 \over 30}A_{4}=\left({7 \over 30}\right)\left({-1 \over 8}\right)A_{0}={-7 \over 240}A_{0}\\[8pt]A_{7}&={3 \over 14}A_{5}=\left({3 \over 14}\right)\left({1 \over 24}\right)A_{1}={1 \over 112}A_{1}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/99e9e4d0385cb69bb0b3738ab539ec825d6911c7)
і ряд розв'язку такий
![{\displaystyle {\begin{aligned}f&=A_{0}z^{0}+A_{1}z^{1}+A_{2}z^{2}+A_{3}z^{3}+A_{4}z^{4}+A_{5}z^{5}+A_{6}z^{6}+A_{7}z^{7}+\cdots \\[8pt]&=A_{0}z^{0}+A_{1}z^{1}+{-1 \over 2}A_{0}z^{2}+{1 \over 6}A_{1}z^{3}+{-1 \over 8}A_{0}z^{4}+{1 \over 24}A_{1}z^{5}+{-7 \over 240}A_{0}z^{6}+{1 \over 112}A_{1}z^{7}+\cdots \\[8pt]&=A_{0}z^{0}+{-1 \over 2}A_{0}z^{2}+{-1 \over 8}A_{0}z^{4}+{-7 \over 240}A_{0}z^{6}+A_{1}z+{1 \over 6}A_{1}z^{3}+{1 \over 24}A_{1}z^{5}+{1 \over 112}A_{1}z^{7}+\cdots \end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2506e46b9672326f09f775182a3e00b93ce99768)
ми можемо розбити у суму двох лінійно незалежних рядів:
![{\displaystyle f=A_{0}\left(1+{-1 \over 2}z^{2}+{-1 \over 8}z^{4}+{-7 \over 240}z^{6}+\cdots \right)+A_{1}\left(z+{1 \over 6}z^{3}+{1 \over 24}z^{5}+{1 \over 112}z^{7}+\cdots \right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3698b9eaaec03138959303391dd4aedcb124c13f)
які можна спростити за допомогою гіпергеометричного ряду.